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Errors in the double variation tecti- 
nique of refractive index measurement 
are analyzed using a new approach. 
The ability to measure matching wave- 
length is characterized, along with the 
effect on the calculated refractive in- 
dex. Refractive index accuracy and pre- 
cision are very dependent on the 
specifics of each calibration set, partic- 
ularly the difference in dispersion be- 
tween the liquid and solid. Our best 
precision (±1 or 2x10"'') is attained 
only when the difference in dispersion 
between liquid and solid is small, and 
is dependent on an individual opera- 
tor's ability to perceive changes in re- 
lief. This precision is impossible to 
achieve for the other glass/liquid combi- 
nations, where we are limited by a pre- 
cision of approximately 1 nm in the 



selection of matching wavelength. A 
bias in the measurement of matching 
wavelength exists that affects the accu- 
racy of the calculated refractive indices. 
The magnitude of the bias appears to 
be controlled by the bandpass of the 
graded interference filter. The errors in 
refractive index using a graded interfer- 
ence filter with a bandpass of 30 nm 
FWHM (full width at half maximum in- 
tensity) are an order of magnitude 
larger than the errors using a filter 
with a bandpass of 15 nm FWHM. 
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1. Introduction 

One of the taslcs in our laboratory was to char- 
acterize the optical properties of asbestos minerals 
to serve as NIST Standard Reference Materials 
(SRMs). Mine-grade chrysotile, amosite, and cro- 
cidolite were to be made available as SRM 1866 [1] 
with certified values for all optical properties, and 
in particular, refractive index. The SRM is in- 
tended to serve as a calibration standard for labo- 
ratories that analyze bulk materials for asbestos 
using polarized light microscopy. Refractive index 
is a primary optical property used to characterize 
transparent minerals, including asbestos, during 
microscopic analysis. Immersion techniques for mi- 
croscopic measurement of refractive index such as 
Becke line, oblique illumination, and focal masking 



are used routinely for refractive index measure- 
ment to the third place [2]. To characterize refer- 
ence materials for use with these techniques, we 
need to use a method with higher accuracy and 
precision. 

There are several microscopic techniques to 
measure refractive index accurately to the fourth 
place, including interferometry [3], the variation of 
temperature at constant wavelength (single varia- 
tion) [4], and the double variation technique [5-8]. 
We decided to use the double variation technique 
for our measurements because we expected it to be 
an improvement on the single variation technique. 
The double variation technique requires the con- 
trol of temperature and wavelength to match the 
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refractive Index of the liquid to that of the un- 
known solid. The solid is immersed in a liquid of 
known refractive index and dispersion, the temper- 
ature is held constant at some value between 20 
and 35 °C, and the wavelength of the illuminating 
light is varied until the refractive index of the solid 
and liquid are observed to match, as indicated by a 
minimum of contrast. The temperature is then set 
to a new value and the process of determining the 
wavelength at which solid and liquid match in re- 
fractive index is repeated. The refractive index of 
the solid is calculated for each set of temperature/ 
matching wavelength measurements using the dis- 
persion equation and temperature coefficient of 
the liquid. The refractive indices thus determined 
are fit to standard dispersion equations, such as the 
Cauchy or Sellmeier equation [9] to describe the 
dispersion of the solid throughout the measured 
range of wavelengths. 

References [6,7] describe a bias in refractive in- 
dex measurements with respect to wavelength us- 
ing the double variation technique which the 
authors attribute to the color sensitivity of the hu- 
man eye. The measured values are biased high at 
wavelengths <555 nm and are biased low at wave- 
lengths >555 nm. The systematic errors reported 
are on the order of ±2x10"". Reference [8] de- 
scribes an approach to correcting the bias by mea- 
suring a glass calibration standard closely matched 
in refractive index to the unknown. A refractive in- 
dex correction value for each wavelength (An^) is 
calculated as nmeas— "true for the glass in question 
and is then applied to the measurements of the 
unknown, resulting in an accuracy and precision of 
approximately ± 1 x lO"*. 

We obtained calibration glasses with refractive 
indices closely matched to the asbestos minerals 
and placed the appropriate glass in the heating 
stage alongside the asbestos mineral as an internal 
standard. We observed a bias with respect to wave- 
length of the same general nature reported earlier, 
however, the magnitude of the bias was at least a 
factor of ten larger. The errors were too large to 
correct for in the manner described in Ref. [8]. We 
began a systematic study of the variables involved 
in the measurement process and developed a dif- 
ferent approach to assessing the errors in refractive 
index. This approach, in which we characterize the 
errors in the determination of the wavelength at 
which solid and liquid match in refractive index at 
each temperature, as opposed to characterizing the 
errors in refractive index with wavelength, provides 
a better understanding of the variables which con- 
trol the accuracy and precision of the technique. 



This error analysis also allowed us to determine 
information concerning the measurement bias and 
the variables that control it. 



2. Calibration 

The calibration procedures for the equipment 
used in the measurement process and the results of 
the calibrations are given below. For a more de- 
tailed description of the types of equipment used in 
the double variation technique, including heating 
stages, illumination sources, immersion liquids, and 
refractometers, see Ref. [6]. 

2.1 Filter Calibration 

We used three different filters in conjunction 
with the quartz halogen light source on our micro- 
scope: 1) a narrow-bandpass graded interference 
filter (GIF), 2) a broad-bandpass GIF, and 3) a set 
of seven fixed-wavelength interference filters. The 
narrow-bandpass GIF has a minimum bandpass of 
15 nm full width at half the maximum intensity 
(FWHM), the broad-bandpass GIF has a minimum 
bandpass of 30 nm FWHM, and the fked- 
wavelength filters each have a bandpass of 10 nm 
FWHM. A GIF is a 20 cm long rectangular inter- 
ference filter that grades in thickness from one end 
to the other, allowing the selection of peak wave- 
lengths from approximately 350 to 750 nm. The 
filter is marked from 1 to 200 in millimeter incre- 
ments, and is calibrated at 16 positions by the man- 
ufacturer to determine the correspondence 
between filter position and peak wavelength. The 
filter holder has an exit slit which can be opened to 
a maximum width of 20 mm. The bandpass of the 
transmission peak is increased by opening the exit 
slit. The set of seven fixed-wavelength filters is de- 
signed to isolate the common spectral lines in the 
visible spectrum with transmittance peaks at 405, 
436, 486, 546, 577, 589, and 656 nm. 

The filters were each calibrated for wavelength 
and bandpass using a spectrophotometer. The 
spectrophotometer was first calibrated for wave- 
length at 589.3 (the mean of the doublet at 589.6 
and 589.0 nm) using a sodium arc lamp, and at 
546.1 and 435.8 nm using a mercury arc lamp. The 
filters were then calibrated using a tungsten light 
source to transmit white light through the filter, 
collecting the transmitted light with the spec- 
trophotometer. The data were corrected for both 
the source characteristics and the relative sensitiv- 
ity of the detector using a blank spectrum collected 
under the same conditions. The peak wavelength of 
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the transmission curve is determined as the cen- 
troid of the peak, and the bandpass is measured 
graphically. The GIFs were calibrated at eight posi- 
tions of the filter, corresponding to a range in 
wavelength of 440-630 nm, with the exit slit fixed at 
2 mm. The narrow-bandpass GIF was also mea- 
sured at one position (wavelength) but with a varia- 
tion in slit width of 2-20 mm. 

Our wavelength measurements for one of the 
GIFs disagreed with the manufacturer's calibration 
measurements by 8 nm, with the same 8 nm dis- 
crepancy at each measured position. The measure- 
ments of the other GIF and the fixed-wavelength 
filters agreed with the manufacturers' values. The 
band passes at FWHM for the two GIFs deter- 
mined with the exit slit at 2 mm agree with the 
manufacturer's values. (The manufacturers do not 
specify a slit size for their bandpass measure- 
ments.) If the exit slit is much narrower than 2 mm, 
the amount of light transmitted is insufficient to 
illuminate the microscope field. The transmission of 
a filter that is graded in thickness is complex, and 
the changes due to changing the size of the exit slit 
are difficult to predict. We measured the transmis- 
sion of the narrow-bandpass GIF as a function of 
slit size, opening the slit symmetrically about the 
center position. The FWHM of the transmission 
peak increased from 15 nm for a 2 mm slit, to 32 nm 
at the fully open position (20 mm). The shape of the 
transmission curve also changed with slit size, with 
the FWHM increasing at a greater rate than the 
width at 1% of the peak height. Therefore, although 
the FWHM of the two GIFs are both approximately 
equivalent when the exit slit of the narrow-bandpass 
GIF is fully opened and the exit slit of the broad- 
bandpass GIF is at the minimum operating width, 
the transmission peaks are not equivalent. 

The two GIFs differ in bandpass (15 or 30 nm 
FWHM) and also in transmission efficiency. The 
manufacturer states a transmission efficiency of 
60% for the narrow-bandpass GIF and 30% for the 
broad-bandpass GIF. We did not measure the 
transmission efficiency of the filters; however, the 
change in brightness of the field of view supports 
the relative difference in efficiency given for the two 
GIFs. We could also see a difference in field illumi- 
nation using the two GIFs that may reflect the dif- 
ference in bandpass. Using the broad-bandpass 
GIF, we could see faint Becke lines with different 
colors than the field, such as red and blue Becke 
lines on a green field, whereas with the narrow- 
bandpass GIF and the fixed-wavelength filters, the 
Becke lines are simply intensity variations of the 
field color. The transmission efficiency of each 



fixed-wavelength filter is listed as approximately 
60%, which again is in relative agreement with our 
observations of brightness in the field of view. 

2.2 Immersion Liquids 

The immersion liquids used in this study were ob- 
tained froni Cargille.' The refractive index at 589.3 
nm was measured on a Zeiss Abbe-type refractome- 
ter (with an Amici prism) and a Bellingham and 
Stanley precision refractometer (without an Amici 
prism). The temperature was measured on the 
Zeiss refractometer with the instrument's fixed 
thermometer, and the temperature on the Belling- 
ham and Stanley refractometer was measured with 
an E-type thermocouple placed on the measuring 
prism. The refractive index measurements of the 
liquids are within ±2x10"'' of the value given by 
Cargille. In addition, the value np-nc, which is a 
measure of the dispersion between 486.1 nm (F) 
and 656.3 nm (C), was determined for each liquid 
on the Zeiss refractometer using the Amici prism, 
and was again found to be within ±2x 10"* of the 
value given by Cargille. Therefore, the dispersion 
equations of the liquids (Cauchy equations) pro- 
vided by Cargille were used in this study. 

We found that the errors in measurements at 
wavelengths other than 589.3 using the Bellingham 
and Stanley refractometer precluded the use of the 
data to determine the dispersion of the liquids inde- 
pendently. The errors were determined by measur- 
ing three SRMs designed for testing refractometers; 
SRM 1822 [10] and SRM 1823 [11]. The question of 
uncertainty in the liquid dispersions is addressed 
later in this paper. 

2.3 Calibration Glasses 

The glasses used in this study were obtained from 
D. Blackburn and D. Kauffman of the Mechanical 
Properties Group at NIST. The glasses were made 
at NIST approximately 30 years ago, and the refrac- 
tive indices were characterized at that time using 
the minimum deviation technique [9] which is a 
high accuracy technique commonly used for the 
measurement of glasses. The technique requires 



' Certain commercial equipment, instruments, or materials are 
identified in this paper to specify adequately the experimental 
procedure. Such identification does not imply recommendation 
or endorsement by the National Institute of Standards and Tech- 
nology, nor does it imply that the materials or equipment identi- 
fied are necessarily the best available for the purpose. 



695 



Volume 97, Number 6, November-December 1992 

Journal of Research of the National Institute of Standards and Technology 



that the glass be in the form of a prism with pol- 
ished surfaces. The nomenclature used to identify 
the glasses when they were prepared was retained 
in this study. Each of the four glasses is identified 
by an alphabetic character followed by three or 
four digits. The refractive indices of some of the 
glasses were measured prior to this study by 
M. Dodge at NIST using the minimum deviation 
technique. The results agree with the earlier mea- 
surements and are accurate to ±5x10"'. The 
values at four wavelengths distributed across the 
visible spectrum were fit to a Cauchy equation us- 
ing a least squares fit. A fit to a linearized 
Sellmeier equation [7] did not differ from a fit to a 
Cauchy equation by more than ±5 x 10"*, which is 
within the level of accuracy of the data; thus the 
mathematically simpler Cauchy equation and fit 
were used. 

2.4 Heating Stage 

The performance of the thermocouple in the 
heating stage during a temperature ramp was de- 
termined by placing a reference thermocouple in 
the immersion cell with its tip in close proximity to 
the stage thermocouple. We detected a problem in 
one of our commercial stages using this procedure. 
We found that the temperature readings were nor- 
mal at room temperature and became anomalously 
high with increasing temperature. The thermocou- 
ple was not out of calibration, but instead was re- 
ceiving additional heat from the wire used to heat 
the metal stage, and therefore was not providing a 
reliable measurement of the temperature of the 
immersion liquid. This problem was rectified by 
changing the placement of the stage thermocouple 
to remove it from proximity to the heating wire. A 
new heating stage was constructed in-house with 
this new design. 

The thermocouple in the heating stage is accu- 
rate to ±0.1 °C at all temperatures between room 
temperature and 35 °C, and there is no tempera- 
ture gradient in the immersion cell. The composi- 
tional stability of the liquids with temperature was 
tested by cycling the temperature up and down 
while measuring the refractive index of the calibra- 
tion glasses. The liquids were found to be stable 
through at least two repeated ramps. 

3. Experimental Design 

We tested the measurement accuracy of the dou- 
ble variation technique using glasses with well 



characterized refractive indices that cover the 
range of the asbestos minerals. There are seven 
sets of glass/liquid calibration data, with a set de- 
fined as the measurements of one glass in one liq- 
uid. The liquids and calibration glasses are listed in 
Table 1 along with the wavelength range over 
which the calibration measurements were per- 
formed. 



Table 1. Calibration sets 



Liquid 


Glass 


Glass 


Glass V 


Am Range 


Id 




«D 




(nm) 


1.552 


F1152 


1.5549 


0.017 


463-564 


1.556 


F1152 


1.5549 


0.017 


526-666 


1.558 


F1152 


1.5549 


0.017 


560-674 


1.674 


A574 


1.6820 


0.018 


497-540 


1.678 


El 889 


1.6783 


0.033 


518-619 


1.682 


A574 


1.6820 


0.018 


545-609 


1.694 


E1442 


1.7110 


0.020 


460-491 



"Dispersive power K=nF-1c/nD-l (Jenkins and White, 
1976). 

Refractive index is measured by placing grains of 
the appropriate glass in the calibrated liquid to se- 
lect the wavelength at which solid and liquid match 
at the given temperature. A match is indicated by a 
minimum of relief. At matching conditions, the dis- 
persion curve of the glass [Eq. (1)] and the disper- 
sion curve of the liquid [Eq. (2)] intersect. Cauchy 
equations, as given in Eqs. (1) and (2), are com- 
monly used to describe the dispersion of materials 
in the visible region of the spectrum [9]. 



rigis-a +"p+"p> 



«iiq=<i+- 



^+4j+[(T- 



[(T-25)-dn/dTl 



(1) 



(2) 



Temperature coefficients (d/i/dT) of the liquids 
are negative, and are on the order of 5 x 10~''/°C. 
Temperature coefficients of most glasses are two 
orders of magnitude lower and are usually 
neglected for this range of temperature. The effect 
of glass temperature coefficients is discussed later 
in the paper. 

For any given intersection (m) of the two disper- 
sion curves, refractive index («), wavelength (A), 
and temperature (T) are uniquely defined. To de- 
termine the accuracy and precision of our measure- 
ments, we need to know the set of (n, A, T)^ for 
r= 20-35 °C, which is our temperature range, for 
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each calibration set. All possible values of rim are 
given by the dispersion equation of the glass. Set- 
ting Eqs. (1) and (2) equal to each other and 
rewriting for T allows Tm to be calculated, as given 
in Eq. (3): 

for rtgis=niiq: 

a+Y2+f,+(25-dn/dT)~(d+f,+^) 



T^ = 



dndT 



(3) 



The resulting (T,\)m data calculated from Eq. (3) 
for each calibration set were fit by least squares to 
a polynomial using DATAPLOT [12] to generate 
an expression that solves for Am at any Tm. For all 
calibration sets, a quadratic or cubic equation was 
sufficient to achieve residuals from the fit of Am to 
Tm of ^ 0.05 nm. 

Because of the relationship among the variables 
rt , A, and T, we can look at our measurement errors 
in terms of any of the three. We chose to analyze 
the errors in T and A, since they are the variables 
that we measure directly, and apply the results to 
the calculation of n. Due to the design of the exper- 
iment, in which we hold temperature constant dur- 
ing the measurement and vary wavelength to make 
a match, it is appropriate to view temperature as 
the independent variable and wavelength as the de- 
pendent variable, and determine our errors in mea- 
suring Am. This assignment of variables fits the 
assumptions of most statistical approaches in which 
there is little or no error in the independent vari- 
able, and all error in the dependent variable. The 
same assumption of error does not hold however, if 
the data are analyzed in the more conventional way 
as (A, n) pairs. 

Figure la illustrates the double variation tech- 
nique and the relationship of n. A, and T in the 
measurement of calibration glass E1442 in liquid 
1.694 from 20 to 35 °C. The intersections between 
the two curves for that range of temperatures are 
shown by the bold line, which contains all possible 
(n, A, T)m for that calibration set. Figure lb is a 
representation of (n. A, 7)^ in A, T space. The 
"true" Am at each T is given by the solid line in 
Fig. lb, with our measurements of A™ given by the 
squares. The Am measurement errors (A Am), calcu- 
lated as AS'^' - Am"°, are shown by the squares in Fig. 
Ic. A linear fit to the measurement errors is shown 
by the solid line in Fig. Ic, and will be discussed 
later in this paper. 

The relationship between the error in the mea- 
surement of Am and errors in the calculated n is 
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Fig. 1. Generation of AAn,: (a) segment of E1442 dispersion 
eurve measured in liquid 1.694 between 20 and 35 °C given by 
bold line, (b) ealibration measurements (squares) vs true (solid 
line), and (e) errors calculated as Ameas-A,r„e (squares) with 
least squares linear fit (solid line). 
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dependent on the difference in dispersion between 
the liquid and the solid. The relationship stems 
from the fact that the refractive index of the solid is 
calculated from the dispersion equation of the liq- 
uid. Figure 2 illustrates this point for a glass mea- 
sured in two liquids of different dispersion. For an 
error of ±1 nm in the determination of Am, the 
error in the calculated n is larger when the glass is 
measured in the high dispersion liquid. This rela- 
tionship between the error in n and the difference 
in dispersion between the glass and the liquid is 
important to the discussion of the technique, and 
will be referred to frequently. We define the differ- 
ence in dispersion between the glass and liquid as 
the difference in refractive index between the glass 
and the liquid at ±1 nm of Am [(«gis-«iiq) @Am±l 
nm]. This value varies slightly with wavelength for 
each glass/liquid calibration set and is a quantita- 
tive description of the relief of the glass grain at 
± 1 nm of match conditions. 



-High dispersion liquid (Lh) 



^Low dispersion liquid (L|) 




Fig. 2. Effect of liquid dispersion on refractive index errors. 

Other variables in the experiment in addition to 
A, n , r, and dispersion include the operator and the 
filter. As mentioned before, the seven calibration 
sets, comprising seven immersion liquids and four 
glasses, were chosen for the analysis of the asbestos 
reference materials, but also allow us to test the ef- 
fects of differences in dispersion. Three operators 
performed the measurements, and differences due 
to operator bias are discussed. The data for the 
seven calibration sets were collected using the nar- 
row-bandpass GIF with a 2 mm exit slit. Additional 
measurements were performed using the fixed- 
wavelength filters, the narrow-bandpass GIF with 
the slit fully opened, and the broad-bandpass GIF 
to test the effects of bandpass. 



Each operator measured at least two grains of 
each glass in each liquid from independent prepara- 
tions. The measurements of Am were performed for 
each glass grain at three temperatures; one at room 
temperature, a second between room temperature 
and 35 °C, and a third at approximately 35 °C. The 
measurements were always made in order of in- 
creasing temperature, and three measurements 
were n\ade at each temperature. The measure- 
ments were made when the temperature in the im- 
mersion cell had stabilized; the temperature did not 
vary by more than 0.1 "C for any set of three mea- 
surements. 

4. Results 

The {T, Am) data collected for each glass/liquid 
calibration set listed in Table 1 were converted to 
(AAm, T), as shown in Figs, lb and Ic, to provide 
common ground on which to compare all calibra- 
tion sets. Both the precision and the accuracy of the 
measurements were determined using the data in 
this form. The accuracy was additionally deter- 
mined by converting the (T, Am) data to (A, n) to 
determine errors in the dispersion curves calculated 
from the measured data. 

4.1 Precision 

The precision of the Am measurement, as indi- 
cated by the variation of the measurements at each 
Tin Figs, lb and Ic, was calculated as one standard 
deviation (1 o-) of the mean of the residuals from a 
linear fit to (AAm, T) for each calibration set. The 
precision was calculated separately for each opera- 
tor's data, and also for the combined set of data. 
The precision of each operator, and of the com- 
bined dataset (all), is given in Table 2, where the 
glass/liquid calibration sets are identified by nu of 
the liquid used. The difference in dispersion be- 
tween liquid and glass is given by («gis -«iiq) @Am ± 1 
nm, which describes quantitatively the relief of the 
glass in the liquid. The slopes of the two curves 
change with wavelength, and therefore (ttgis-nnq) 
@An,±l nm varies with wavelength. The range 
given for (ngis -niiq) @Am ± 1 nm for each calibration 
set represents a dependence on wavelength, and de- 
creases with increasing wavelength. The range in Am 
for each calibration set is given in Table 1. 

The precision can be interpreted with respect to 
Fig. 2, in which the solid is measured in either a 
high dispersion liquid (V> 0.050) or a low disper- 
sion liquid (F = 0.031). (The dispersion of the glass 
remains relatively constant (K = 0.017 -0.020), 
except for calibration set 1.678 which has a higher 
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Table 2. 


A„ precision 










Set 


Dpi 


Op2 

1(tA„ 


Op3 

(nm) 


All" 


LiqK" 


("gls-«liq) 

@A„±1 nm 
xlO-" 


1.552 


±2.8 


±3.7 


±2.2 


±3.3 


0.031 


0.9-0.8 


1.556 


±2.8 


±3.5 


±3.1 


±4.6 


0.031 


0.6-0.3 


1.558 


±3.1 


±3.5 


±3.2 


±3.5 


0.031 


0.5-0.3 


1.674 


±1.7 


±2.0 


±1.3 


±2.0 


0.050 


1.9-1.7 


1.678 


±1.5 


±2.9 


±1.1 


±2.0 


0.050 


0.9-0.6 


1.682 


±1.6 


±2.7 


±1.7 


±2.3 


0.053 


1.4-1.0 


1.694 


±1.6 


±2.9 


±1.3 


±2.4 


0.053 


3.1-2.3 



' Determined from eombined dataset. 
'' Dispersive power K=«f-'Jc/«d-1. 



dispersion glass (F = 0.033).) The Am measure- 
ments are more precise in the high dispersion liq- 
uids than in the low dispersion liquids. The 
precision is operator dependent, but this basic cat- 
egorization holds for all three. Operators 1 and 3 
have an average lo-of approximately 1.5 nm for the 
high dispersion liquids and an average lo-of ap- 
proximately 3 nm for the low dispersion liquids. 
Operator 2 is less precise but follows the same 
trends. The improvement in precision with increas- 
ing liquid dispersion and therefore increasing 
(«gis-niiq) @Am± 1 nm indicates that the amount of 
relief influences the precision of the measurement, 
which is a reasonable conclusion. The precision 
also seems to have a lower boundary at la- equal to 
approximately 1 nm, beyond which the relief can 
increase without an additional improvement in pre- 
cision. The precision for set 1.694 is comparable to 
the other high dispersion liquid sets, even though 
this set has the highest relief. This probably repre- 
sents a limitation of the GIF in that it may not be 
capable of providing peak wavelengths with a sepa- 
ration of better than approximately 1 nm. 

The precision in n for the measurements can be 
calculated by multiplying the A™ precision by 
(«gis - rtiiq) @Am ± 1 nm, which defines the error in n 
for each 1 nm error in A^. The precision in n , cal- 
culated using the mean value of («gis-«iiq) @Am± 1 
nm for each calibration set, is given in Table 3. The 
precision in n for the low dispersion liquids, where 
relief is the controlling factor, is be'tween 
±1x10"" and ±2.4x10'" for operators 1 and 3. 
The precision degrades for the high dispersion liq- 
uids with the exception of set 1.678, which has a 
high glass dispersion and thus a lower relief. The 
precision in Am, and thus n , for the high relief cali- 
bration sets is ultimately controlled by the ~1 nm 
limitation in Am. Therefore, for calibration set 
1.694, the precision in n is poor even though the 
precision in Am is comparable to other sets, due 



solely to the high relief. Put in another way, to 
achieve a precision in n of approximately 
± 1 X 10"" for calibration set 1.694 requires a preci- 
sion in the measurement of Am of approximately 
± 0.4 nm or better, which appears to be beyond our 
measurement capabilities. 

The precision of the measurements can be sum- 
marized as follows: at best we can discriminate 
changes in relief of approximately 1 X 10"" and our 
best precision in the measurement of Am is ~ 1 nm. 
Because we cannot measure A™ with a precision 
better than ~ 1 nm, a difference in dispersion be- 
tween the liquid and solid which results in 
(«gis-«iiq) @Am±l nm of greater than ±1x10"", 
will result in an imprecision in n that is larger than 
±1x10"". 

The precision of the Am measurements, as given 
in Table 2, is described in terms of individual oper- 
ators, and also in terms of the combined data set. 
For each calibration set, lo-for the combined data- 
set is larger than the average of the three opera- 
tors, indicating a bias among operators. Analysis of 
the average value of Am at a given temperature for 
each operator indicates that operator 1 has a con- 
sistent negative bias on the order of approximately 
1 nm with respect to the other two operators. It is 
possible to correct this type of problem through re- 
training of the operator or by the determination of 
calibration curves for each operator. 

Table 3. n measurement preeision 



Set 


Dpi 


Op 2 
\(rn (xlO-") 


Op 3 


All 


1.552 


±2.4 


±3.1 


±1.9 


±2.8 


1.556 


±1.3 


±1.6 


±1.4 


±2.1 


1.558 


±1.2 


±1.4 


±1.3 


±1.4 


1.674 


±3.1 


±3.6 


±2.3 


±3.6 


1.678 


±1.1 


±2.2 


±0.8 


±1.2 


1.682 


±1.9 


±3.2 


±2.0 


±2.8 


1.694 


±4.3 


±7.8 


±3.5 


±6.5 


4.2 


Accuracy 









The accuracy of the technique was analyzed by 
determining the errors in the measurement of Am 
with respect to T for each glass/liquid calibration 
set. Systematic errors in the measurements of Am 
were observed for each calibration set. The linear 
fits of AAm to T for each set have negative slopes, 
and the values of AAm predicted by the fits are pos- 
itive or close to zero at the low temperature end of 
each set, and negative at the high temperature end 
of each set, as shown in Fig. 3a. Because of the 
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nature of the experiment, an increase in tempera- 
ture always corresponds to a decrease in Am, and 
therefore a bias associated with temperature can- 
not be separated from a bias associated with wave- 
length. Linear fits to the errors in Am with respect 
to wavelength for each set are given in Fig. 3b, with 
temperature listed at each end point. The linear 
fits all have positive slopes, and the predicted 
AAm's are negative at the short wavelength end of 
each set, and positive at the long wavelength end of 
each set. It is important to note that the calibration 
sets cover different ranges of wavelength, and that 
the errors are not associated with wavelength in an 
absolute sense. These data indicate that the sys- 
tematic errors are not due to color sensitivity of the 
detector (human), as we will discuss later in the 
paper. The linear fits to A Am with respect to tem- 
perature and wavelength shown in Fig. 3 are noisy, 
as can be seen in Fig. Ic, and the absolute values of 
the slope and intercept for each set are not the 
same. However, the repetition of the trends in the 
errors for all sets is significant, and indicates that 
there is a bias in our measurements which corre- 
lates with temperature and wavelength. 

Some possible sources of error that could pro- 
duce the observed bias include those related to the 
heating stage (thermocouple calibration, tempera- 
ture instability, temperature gradients), systematic 
errors in the calibration of the GIF, systematic er- 
rors in the dispersion of the liquids and/or glasses, 
and uncertainties in temperature coefficients. The 
calibrations of the heating stage, GIF, and glasses 
are described in the calibration section and we can 
eliminate these factors as possible significant con- 
tributors to the observed bias. We were concerned 
about possible systematic errors in the calibration 
of the liquids, since we were not able to measure 
the dispersion independently (other than the mea- 
surement of F-C). The dispersions of the liquids 
also change with temperature such that dn/dT is 
not the same at all wavelengths, although by a 
small amount (10~^). In addition, the temperature 
coefficients of the glasses are not known, although 
they are also small (lO"*). 

The possibility that systematic errors in the 
dispersions of the liquids and uncertainties in tem- 
perature coefficients were responsible for the bias 
observed in Am was tested by using a liquid and 
glass for which those quantities are well character- 
ized. We measured SRM 1822, an optical glass re- 
fractive index standard, in SRM 18231, one of the 
refractive index liquid standards. The dispersions 
of both liquid and glass are well characterized, as 
are the temperature coefficients of the liquid with 
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Fig. 3. Bias in A„; (a) with temperature and (b) with wave- 
length. Calibration sets identified by no of liquid. Lines repre- 
sent least squares linear fits to A A™ for each calibration set. 

respect to wavelength. SRM 18231 is a silicone oil 
that is chemically and thermally stable. The refrac- 
tive index of the liquid was measured using the 
minimum deviation technique and is certified at 10 
wavelengths in the visible with an accuracy of 
±4x10"^ at four temperatures from 20 to 80 °C. 
SRM 1822 is a commercial soda-lime glass that is 
certified for refractive index at 13 wavelengths in 
the visible with an accuracy of ±9x10"*. The 
temperature coefficient for SRM 1822 is not given 
on the certificate, however, the temperature 
coefficients for common optical glasses are avail- 
able from the Schott catalog [13] and range from 
- 3 X lO"" to - 5 X 10"*, increasing with increasing 
refractive index. A barium crown glass (BK7) made 
by Schott Glass is similar in refractive index to 
SRM 1822, and has a temperature coefficient of 
-3x10"*. This value was taken as the tempera- 
ture coefficient for SRM 1822, although the 
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selection of other temperature coefficients within 
the range listed produced no significant difference 
in results. 

The errors in the measurements of SRM 1822 in 
SRM 18231 using the narrow bandpass GIF are 
shown in Fig. 4, and the linear fit to the errors 
follows the same trend as the fits to the errors for 
all other calibration sets. The conclusion from 
these measurements is that the bias in the mea- 
surements of our seven calibration sets, as dis- 
played in Fig. 3, is not due to uncertainties in dis- 
persion, temperature coefficients, or liquid 
instability at high temperatures. 

We were able to change the magnitude of the 
bias observed in Am by changing the bandpass of 
the GIF. We performed measurements of SRM 
1822 using the narrow bandpass GIF with the slit 
fully open (FWHM —30 nm), and the broad band- 
pass GIF, and the errors are illustrated in Fig. 4. 
The linear fits to the data from the GIFs have neg- 
ative slopes, and the magnitude of the slope in- 
creases with increasing bandpass. The slopes of the 
three linear fits are significantly different from 
each other such that at 30 "C the 95% confidence 
limits for each fit comprise separate populations. 
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Fig. 4. Effects of bandpass on measurement bias for measure- 
ments of SRM 1822 in SRM 18231. The broad-bandpass GIF 
(squares) has a 30 nm FWHM, the narrow-bandpass GIF has a 
15 nm FWHM (stars) which increases to 30 nm with the slit fully 
open (triangles). The fixed-wavelength filters each have a 10 nm 
FWHM (diamonds). 

We do not see a bias with temperature or wave- 
length when we use the fixed wavelength filters 
(10 nm FWHM) and vary the temperature to mea- 
sure Tm (single variation technique). Measurements 
of SRM 1822 using the 546.1, 486.1, and 435.8 nm 
filters are shown in Fig. 4 [after conversion to (7, 
A Am)] along with a least squares linear fit. A 95% 



confidence interval of the slope of the fit includes 
zero, indicating that there is no significant bias for 
this data. 

Measurements of the other calibration sets using 
the fixed-wavelength filters also do not exhibit any 
systematic errors. Due to the spacing of the trans- 
mittance peaks of the filters, we are usually limited 
to one or two temperature measurements for each 
calibration set, and therefore, we have smaller data 
sets than for the GIFs. However, the measurements 
from the fixed-wavelength filters span the full tem- 
perature and wavelength range, and we do not ob- 
serve a bias with respect to temperature or wave- 
length. This result conflicts with that of Ref. [7] 
which reports the same refractive index bias for 
both fixed-wavelength filters and a GIF. The au- 
thors do not state the bandpass of their filters, and 
it is possible that the different results are due to this 
factor. We did not test any fixed-wavelength filters 
with a broader bandpass. We did test the possibility 
that the variation of wavelength was responsible for 
the bias by fixing the position of the GIF and per- 
forming a single variation measurement, as with the 
fixed-wavelength filters. The data collected with the 
GIF in a fked position were biased in the same 
fashion as the data collected by varying the wave- 
length. 

4.3 Dispersion Calculations 

Given the fact that there is a bias in the measure- 
ment of Am, we need to determine the effect of this 
bias on the calculated refractive indices and fitted 
dispersion curves. Qualitatively, the errors in n 
always have the opposite sign of A Am, therefore An 
is positive for the short wavelength— high tempera- 
ture end of each calibration set, and negative for 
the long wavelength— low temperature end of each 
calibration set. We can get a general idea of the 
magnitude of A/j using the same procedure we used 
for determining the precision, namely by multiply- 
ing the error in Am by (rtgis-niiq)@Am± I nm for that 
calibration set. For example, from Fig. 3b we can 
see that at the short wavelength end of each calibra- 
tion set, the predicted errors in Am from the linear 
fits range from approximately -2 nm to -6 nm, 
and we will assume an average error of - 4 nm for 
discussion purposes only. If we multiply — 4 nm by 
the largest value of (ngi5-niiq)@Am± 1 nm for each 
set from Table 2 (which corresponds to the short 
wavelength end) we see that the errors in n will 
range from + 2 x 10"" to -I- 1.2 X 10"^ As discussed 
with reference to precision, the errors in the 
calculated n increase as the difference in disper- 
sion between the glass and the liquid increase, even 
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for the same error in Am. This fundamental associa- 
tion stems from the assumption of the technique 
that the measured value lies on the liquid disper- 
sion curve. 

The ultimate goal of the measurement technique 
is to calculate dispersion constants for the un- 
known to predict the refractive index throughout 
the measurement range. We must, therefore, de- 
termine the errors involved in the calculation of 
the dispersion of the calibration glasses. The dis- 
persion of each glass listed in Table 1 was calcu- 
lated by converting the measurements from each 
liquid to (A, «), combining the data from each liq- 
uid, and fitting the combined set to a Cauchy equa- 
tion. In the case of glass F1152 there are three sets 
of data which are combined for the fit, for glass 
A574 there are two sets of data, and for glasses 
E1889 and E1442 there is only one set of data 
each. The (A, n) fits to the data for each glass are 
shown in Fig. 5, along with the "true" (A, n) shown 
by the bold line. The error in n at each A is deter- 
mined by subtracting the "true" n at each wave- 
length from the n predicted by the fit to the data. 
The Art's thus calculated for the glasses are shown 
in Fig. 6. This is the same approach recommended 
by Su et al. (1987), who use three liquids to mea- 
sure each calibration glass. 

The errors in refractive index calculated using 
this approach now exhibit the systematic bias dis- 
cussed by Refs. [6-8], particularly for glasses F1152 
and A574, which are biased high for short wave- 
lengths and biased low for long wavelengths, ap- 
proaching zero error at approximately 550 nm. As 
discussed earlier, however, the bias does not corre- 
late with absolute wavelength, but with the mini- 
mum and maximum wavelength of each glass/liquid 
calibration set. The apparent correlation of the 
bias with absolute wavelength for glasses F1152 
and A574 results from the combination of the data 
from multiple liquids, and the overlap of the sets. If 
the data from each calibration set for glass F1152 
are treated separately, we obtain the Awa curves 
shown by the thin lines in Fig. 7, whereas if we 
combine the sets we obtain the An a curve shown by 
the bold line in Fig. 7. The error in refractive index 
at a particular wavelength is dependent upon the 
number of liquids used to measure the glass; the 
bias in each set of data from a given liquid will be 
compensated over the wavelength range where 
there is overlap with data from another liquid. 

The dependence of An a on the number of liquids 
used in the calibration measurements and on the 
difference in dispersion between the liquid and the 
solid indicates that the use of such values to correct 
measurements of unknowns is only appropriate 



when the unknown has the same dispersion as the 
calibration material (or the difference in dispersion 
between liquid and solid is the same for both) and 
when the same number of liquids covering the 
same general range of wavelengths are used. For 
example, if only one liquid were used to determine 
AnsM for glass F1152, the result would be either 
- 1 X 10""*, 0, or + 1.5 X 10"*, depending on the liq- 
uid used for calibration. In addition, glass E1889 
would be an inappropriate calibration material to 
use to correct the measurements of glass A574 
even though they are similar in refractive index, 
due to the large difference in dispersion. 

The refractive index errors are also highly de- 
pendent on the filter used because of the bias asso- 
ciated with bandpass. Figure 8 shows the calculated 
errors in n for glass SRM 1822, using the data from 
Fig. 4. The error in refractive index for the fixed 
wavelength filters and the narrow bandpass GIF 
with the 2 mm slit is within ±2x10"*, but 
increases to approximately + 1 x 10"^ for the broad 
bandpass GIF. 

5. Discussion 

We believe that the bias we observe in Am may 
explain the bias observed in refractive index in the 
previous studies [6-8]. This bias correlates with 
temperature and wavelength, but cannot be traced 
to calibration or measurement errors associated 
with either variable. The magnitude of the bias ap- 
pears to be directly influenced by the bandpass of 
the interference filter. We do not observe a bias 
with the 10 nm FWHM fixed wavelength filters, but 
we observe a bias with the variable wavelength 
GIFs that increases with increasing bandpass. We 
did not test fixed wavelength filters with wider 
bandpasses and the question remains as to whether 
the GIFs themselves are a source of the bias, with 
the bandpass of the GIFs as an additional factor, 
or whether it is bandpass alone (or another vari- 
able associated with the transmission of the filter) 
which is the critical factor. The GIFs and the fixed 
wavelength filters are similar in that both are inter- 
ference filters, but are dissimilar in transmission 
characteristics, due to the gradation in thickness of 
the GIF. Louisnathan et al. (1978) report a bias 
using fixed wavelength filters of unspecified band- 
pass, and it is therefore probable that the bias is 
not restricted to the use of GIFs. We are perplexed 
by the apparent association between bandpass and 
the observed bias, and can only suggest that the 
relief of the glass grains must be affected by the 
bandpass in some manner that biases our selection 
of A„. 
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Fig. 5. Measurements of calibration glasses with fitted Cauchy (thin line) vs true (bold line). 
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Fig. 6. Dispersion curve measurement errors for calibration 
glasses. 
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Fig. 7. Dispersion curve errors from multiple vs single liquids. 
Dispersion curve measurement errors for glass F1152 deter- 
mined for each liquid (thin lines) and for combined dataset 
(bold line). 
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Fig. 8. Effect of bandpass on refractive index errors using data 
displayed in Fig. 4. 



We found that even with the bias in the tech- 
nique, the measurements of the glasses were ac- 
ceptable (within ±5x10"*) with the exception of 
glass E1442, for which we have a maximum error of 
-f 1.2x10"^ and a precision of ±6x10"*. The 
large error in refractive index is due to the large 
difference in dispersion between the glass and the 
liquid, and not to any degradation in the actual 
measurements. This calibration set emphasizes the 
need to select liquids with the lowest dispersion 
possible although the choices are usually quite lim- 
ited. In the initial stages of this project, when we 
were determining accuracy and precision in the 
conventional manner by analyzing the errors in re- 
fractive index with wavelength, we mistakenly 
thought that the precision and accuracy of the 
technique were dependent on the absolute refrac- 
tive index of the glass as our errors increased with 
increasing refractive index. Only after we adopted 
the approach of assessing the errors in Am did we 
realize that the correlation was with the difference 
in dispersion between the liquid and the glass, 
which generally increases with increasing refractive 
index. 

Because of the general acceptability of the dou- 
ble variation measurements, we felt it was appro- 
priate to use the technique for the measurement of 
the asbestos reference materials. We could have 
used the fixed-wavelength filters and the single 
variation technique, for which we did not detect a 
bias, but our heating stage is not designed for rapid 
oscillation and equilibration of temperature which 
is necessary for efficient single variation measure- 
ments. We do not correct the measurements of our 
unknowns using the calibration measurements be- 
cause of the strong dependence of Anx on the 
specifics of each calibration set. We do use the cal- 
ibration measurements to provide estimates of the 
errors in our measurements and to determine 
whether our variables are under control during 
measurement of the unknowns. 

6. Conclusions 

The errors in refractive index measurement us- 
ing the double variation technique were deter- 
mined by characterizing the accuracy and precision 
of measuring the matching wavelength (Am). The 
precision and accuracy of the technique in terms of 
the measurement of refractive index are ultimately 
dependent on the difference in dispersion between 
the solid and the liquid. The best precision possi- 
ble, which in our case is 1 or 2x 10"'', is dependent 
on the operator's ability to perceive changes in 
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relief, and is only possible for those liquid/solid 
combinations where (ngis-niiq)@Am + l nm < 
1 X 10"''. This limitation is based on our inability to 
measure Am with a precision better than approxi- 
mately ±1 nm, which may represent a limitation of 
the filter. 

There is a bias in our measurements of Am using 
a GIF which correlates with temperature and 
wavelength. The bias results in positive refractive 
index errors at the short wavelength— high temper- 
ature end of each dataset, and negative refractive 
index errors at the long wavelength— low tempera- 
ture end of each dataset. (A dataset is defined as 
the measurements of 1 glass in 1 liquid). The bias 
does not correlate with absolute wavelength, and 
therefore does not appear to result from a relative 
sensitivity of the human eye. The magnitude of the 
bias increases with increasing bandpass, as deter- 
mined by opening the exit slit of the narrow-band- 
pass GIF, and by using a GIF with a broader 
bandpass. The bias, in association with a glass and 
liquid with a large difference in dispersion, can 
result in errors in refractive index that exceed 
±5x10"*. 

In general, our errors in refractive index using 
the narrow-bandpass GIF (15 nm FWHM) are 
within ±5x10"'' throughout the visible spectrum, 
and within ±1x10"* at 589.3 nm, as calculated 
from fits to the combined data from multiple liq- 
uids. Correction factors calculated from the mea- 
surements of calibration glasses are dependent on 
the number of liquids used and the dispersion of 
both liquid and solid, and should ideally be used 
only when these characteristics are matched in the 
measurement of the unknown. 
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